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ABSTRACT. We investigate the nature of instabilities in cosmic ray 
shocks by using two distinct models for the shock wave. For wavelengths 
which are short relative to the thickness of the shock wave, the shock is 
treated as a smoothly decelerating flow, and an appropriate JWKB type ex- 
pansion is used to describe the perturbations to the flow. In this, the^ 
short wavelength regime, the presence of "squeezing" and an effective "g" 
renders strong cosmic ray shocks unstable in a way which is similar to 
instabilities in other supersonic flows, such as in de Laval nozzle flow 
or a heat conduction dominated shock wave. In the long wavelength limit, 
where the shock is treated as a discontinuous transition, we derive a 
"stability function" which, if negative, corresponds to unstable dis- 
turbances growing exponentially in time. In this case, we find that if 
the cosmic ray fluid is relativistic (y c = V and the background plasma 
ideal (y = %), then strong shocks are unstable. 

1. INTRODUCTION. We examine the stability of cosmic ray shocks (Axford 
et al [1977], Drury and VSlk [1981], Axford et al [1982], Vblk [1984], 
Drury [1984], Dorfi [1984]) in both short and long wavelength regimes. 

The short wavelength instability in the two-fluid description of a cosmic 
ray shock may be shown to be similar to the instabilities in other super- 
sonic decelerating flows. Erpenbeck, [1962], established that discon- 
tinuous transitions (shocks) are always stable for ideal fluids. We 
generalize his results, and show that the relativistic two-fluid model 
considered here is unstable. 

2. THE GOVERNING EQUATIONS. The equations which describe the self- 
consistent interaction between a background thermal plasma and a 
relativistic cosmic ray fluid are given by 
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where p, p, p c , and u denote gas density, gas pressure, cosmic ray pres- 
sure and gas velocity respectively, and k an effective cosmic ray 
diffusion coefficient. These equations are supplemented by the cosmic 
ray energy equation 

h ( J^-j) + yFj div tY c P c U " <V Pc ] = U. V p c . 
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which, together with 2.1(a,b), enables 2.1(c) to be expressed in the 
adiabatic form 
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(2.3) 


where ^ = . "ft + ii- v is the convective derivative. The derivation 
?iQ?m SCUSS10n ^ ese equations may be found in, for example, Drury, 


A perturbation analysis of (2.1), (2.2) and (2.3) shows that short 
wavelength perturbations propagate at the plasma sound speed, c, where 



(2.4) 


and hence decouple from the cosmic ray fluid. Conversely, long wave- 
length disturbances propagate at the mixed sound speed. 


Ptsukin, [1981], 



^Vc + V /p » 


a c * cosmic ray fluid sound speed. 


(2.5) 


A discussion of the shock structure problem may be found in Axford et al , 
[1982], and Drury and Volk, [1981], McKenzie and Volk [1982], Volk et al 


3. THE SHORT WAVELENGTH INSTABILITY. For the sake of completeness we 
present a brief modified account of the instability discussed by Drury 
[1984]. By considering waves whose wavelengths are short relative to 
the thickness of the shock wave, one may treat the transition as a 
continuous decelerating flow. If one views the decelerating flow as an 
"effective gravitational force", the stability problem has the character 
of a "generalized Rayleigh-Taylor" type problem for a disturbance 
propagating through the system. The introduction of compressive pertur- 
bations = (p» P> Pc» u) on the background state t[> 0 = (p 0 , p 0 , p co , u 0 ) 
enables equations (2.1) and (2.3) to be linearized. We introduce a shock 
thickness parameter L, related to the decelerating flow u Q , by defining 

1 - -J_ du 0 

L u o dx ' 

, du n 

where g* is the flow velocity profile. On either side of the shock, 
u 0 is constant, and hence X = 0. By using the steady state conservation 
laws (2.1) and (2.3), one L can express the gradients of the remaining 
background variables p 0 , p 0 , p co , as functions of L and the Mach number 
of the flow. We restrict the analysis to velocity profiles which vary 
slowly in space, and thus need only consider terms which are 0(X). This 
is tantamount to performing a WKBJ expansion in space. The use'-of this 
expansion, and the search for plane wave solutions whose amplitudes 
vary harmonically, reduces (2.1) - (2.3) to the following dispersion 
equation; 
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co' + iu 0 
L 

(3.2) 

where to' = co - u 0 k in the Doppler shifted frequency, and k the wave num- 
ber. In the absence of a shock, the dispersion relation reduces to 

co' = 0; to 1 = ±c 0 k. (entropy-vorticity; sound waves) 

A perturbation analysis of (3.2), about these frequencies, indicates 
that in the presence of a decelerating flow, 

Im co' = - f 0 . ( 3 - 3 ) 

for an entropy wave, and for a sound wave 

l m w . = ^ - E£ 2 t M 2 - (y+1) M - 1], (3.4) 

which, since we consider solutions of the form expicot, implies insta- 
bility for sufficiently large Mach numbers M. 

4. THE LONG WAVELENGTH INSTABILITY. In this section we generalize 
and extend the results of D' Iakov, [1958], Erpenbeck, [1962 a,b], and 
McKenzie and Westphal , [1968] to include the effect of a cosmic ray 
fluid. We derive the transmission coefficient of an acoustic wave 
incident on an oblique shock, and examine the singularities of this 
transmission function. In the case of a supersonic-subsonic transition 
it is straightforward to show that, if the shock is perturbed, then 
three waves diverge from the shock (and thus the problem is well posed). 
After perturbing the generalized Rankine-Hugoniot boundary conditions 


[pu n ] = 0; (continuity) (4.1a) 

[p + P c + pu n ] = 0; (normal momentum) (4.1b) 

[u ] = 0; (tangential momentum) (4.1c) 

[|u n 2 + _L_ £ + Yc_ Pp] = 0, (energy) (4. Id) 

Y_i P Yc-1 P 

where n and t represent directions normal and tangential to the shock, 
one may derive the transmission coefficient for an acoustic wave, 

Sp 2 a 3 (1 - ^Al) Ni + MiN 2 + $ 3 N 3 , (4.2) 

6 Pl a 3 (l- 1W) Dj + 8 3 6 2 D 2 - 6 3 D 3 


where (a., 3-) denote the cosine and sine of the angle between Jc and the 
x-axis, Jnd n Ni, D-; are suitable generalizations of the functions presen- 
ted in McKenzie and Westphal, [1968]. The singularities of the 
function (4.2) correspond to the dispersion equation. 
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On defining the variables 



Y-l 

R* = 1 + Vi 


ac 2 1 R 2 

c » M = M g 2 , M $ = gas Mach number, 


the stability function, F , may be written as 


F s ~ p7 0-M 2 )(2 -(y-1) -jr* (p 2 -1)) - (1 + (y-1)!t! ~ ( 1+ (Y-1) 



(4.3) 

for waves propagating with the flow, and where the subscript 1(2) 
indicates ahead (behind) the shock. In the cold plasma approximation, 
cosmic ray shocks are unstable for all supersonic - subsonic transitions. 
In the alternate approximation of a tenuous cosmic ray gas ahead of the 
shock, it may be established that shocks are unstable for all Mach 
numbers £ 1,3. The growth rate of these instabilities (for strong 

shocks) is of the order the compression ratio (£ 2 ) x ku x . 
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